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■ Abstract 

The gauge interactions of any supersymmetric extension of the stan- 
dard model involve new flavor mixing matrices. The assumptions involved 
in the construction of minimal supersymmetric models, both SU(3) x 
SU(2) x U(l) and grand unified theories, force a large degree of triviality 
on these matrices. However, the requirement of realistic quark and lepton 
masses in supersymmetric grand unified theories forces these matrices to 
be non-trivial. This leads to important new dominant contributions to 
the neutron electric dipole moment and to the decay mode p — * K°fi + , 
and suggests that there may be important weak scale radiative correc- 
tions to the Yukawa coupling matrix of the up quarks. The lepton flavor 
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violating signal ji — > is studied in these theories when tan/3 is suf- 
ficiently large that radiative effects of couplings other than X t must be 
included. The naive expectation that large tan/3 will force sleptons to 
unacceptably large masses is not borne out: radiative suppressions to the 
leptonic flavor mixing angles allow regions where the sleptons are as light 
as 300 GeV, provided the top Yukawa coupling in the unified theory is 
near the minimal value consistent with m<. 
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I. Introduction 

It has recently been demonstrated that flavor and CP violation provide 
an important new probe of supersymmetric grand unified theories [1-4]. These 
new signals, such as \i — > ej and the electron electric dipole moment d e , are 
complementary to the classic tests of proton decay, neutrino masses and quark 
and charged lepton mass relations. The classic tests are very dependent on the 
flavor interactions and symmetry breaking sector of the unified model: it is only 
too easy to construct models in which these signals are absent or unobservable. 
However, they are insensitive to the hardness scale, Ah, of supersymmetry break- 
ing.Q On the other hand, the new flavor and CP violating signals are relatively 
insensitive to the form of the flavor interactions and unified gauge symmetry 
breaking, but are absent if the hardness scale, Ah, falls beneath the unified 
scale, Mq- The signals are generated by the unified flavor interactions leaving 
an imprint on the form of the soft supersymmetry breaking operators ||, which 
is only possible if supersymmetry breaking is present in the unified theory at 
scales above Mq. 

The flavor and CP violating signals have been computed in the minimal 
SU(5) and 5*0(10) models for leptonic [1-3] and hadronic processes f§, for 
moderate values of tan/?, the ratio of the two Higgs vacuum expectation values. 
While rare muon decays provide an important probe of SU(5), it is the SO(10) 
theory which is most powerfully tested. If the hardness scale for supersymmetry 
breaking is large enough, as in the popular supergravity models, it may be 
possible for the minimal SO(10) theory to be probed throughout the interesting 
range of superpartner masses by searches for /i — > and d e . 

The flavor changing and CP violating probes of 5*0(10) are sufficiently 
powerful to warrant an exploration of consequences for non-minimal models, 
which is the subject of this paper. In particular, we study 50(10) theories in 
which 

(I) The Yukawa interactions are non-minimal. 

In the minimal model the quarks and leptons lie in three 16 's and the 
two Higgs doublets Hjj and Hp lie in two 10 dimensional representations 10^ 

1 This is the highest scale at which supersymmetry breaking squark and gluino masses 
appear in the theory as local interactions. 



1 



and 10 d- The quark and charged lepton masses are assumed to arise from the 
interactions 16A[/16 10u + 16Aj)16 10d- This model is a useful fiction: it is very 
simple to work with, but leads to the mass relation m e jm ll = m^/ms, which is in 
error by an order of magnitude. It is clearly necessary to introduce a mechanism 
to insert SO(10) breaking into the Yukawa interactions. The simplest way to 
achieve this is to assume that at the unification scale, Mq, some of the Yukawa 
interactions arise from higher dimensional operators involving fields A which 
break the SO (10) symmetry group. This implies that Xu,d —> ^u,d(A). Every 
realistic model of SO (10) which has been constructed has this form; hence one 
should view this generalization of the minimal model as a necessity. 

(II) The ratio of electroweak VEV's, tan/5 = vu/vd, is allowed to be large, 
« nit/nit,. 

This is certainly not a necessity; to the contrary, a simple extrapolation of 
the results of @ to such large values of tan (3 suggests that it is already excluded 
by the present limit on fi — > &y. The case of large tan/3 in SO (10) has received 
much attention |], [7|, |9|] partly because it has important ramifications for the 
origin of m t /mb = (Xt/Xb) tan/3. To what extent is this puzzling large ratio to 
be understood as a large hierarchy of Yukawa couplings, and to what extent in 
terms of a large value for tan/3? If the third generation masses arise from a 
single interaction of the form I63I63IO it is possible to predict m t using and 
m T as input ||, providing the theory is perturbative up to Mq. The prediction 
is 175 ± 10 GeV 0, and requires tan/3 ~ mt/mb- In this paper we investigate 
whether this intriguing possibility is excluded by the fi — ► ej signal; or, more 
correctly, we determine whether it requires a soft origin for supersymmmetry 



breaking, making it incompatible with the standard supergravity scenario |10 



In the next section we show that SO (10) models with A — > X(A) possess 
new gaugino mixing matrices in the up-quark sector, which did not arise in 
the minimal models. In section III we set our notation for the super symmetric 
standard model with arbitrary gaugino mixing matrices, and we show which 
mixing matrices are expected from unified models according to the gauge group 
and the value of tan /3. In section IV we describe the new phenomenological 
signatures which are generated by the gaugino mixing matrices in the up sector; 
these signatures are generic to all models with Yukawa interactions generated 
from higher dimensional operators. The consequences of large tan/3 for the 
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flavor and CP violating signatures are analyzed analytically in section V and 
numerically in section VI. The analysis of the first five sections applies to a wide 
class of models. In section VII we illustrate the results in the particular models 
introduced by Anderson et al |J . As well as providing illustrations, these models 
have features unique to themselves. Conclusions are drawn in section VIII. 
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II. New Flavor Mixing in the Up Sector 

In [1-4] flavor and CP violating signals are studied in minimal SU(5) and 
50(10) models with moderate tan/3. In these models the radiative corrections 
to the scalar mass matrices are dominated by the top quark Yukawa coupling At 
of the unified theory, so the scalar mass matrices tend to align with the up-type 
Yukawa coupling matrix and all non-trivial flavor mixing matrices are simply 
related to the KM matrix. However, as mentioned above, the minimal models 
do not give realistic fermion masses. One has to insert 50(10) breaking into 
the Yukawa interactions. The simplest way to achieve this is to assume that the 
light fermion masses come from the non-renormalizable operators 

Ai Ao At A t+1 A n , . 

A « 16 'MTl|-i 10 A^'"lt 16 - (21) 

where the 16j's contain the three low energy families, 10 contains the Higgs dou- 
blets, and A's are adjoint fields with vacuum expectation values (VEV's) which 
break the 5O(10) gauge group. After substituting in the VEV's of the adjoints, 
they become the usual Yukawa interactions with different Clebsch factors as- 
sociated with Yukawa couplings of fields with different quantum numbers. For 
example in the models introduced by Anderson et al. ||, (hereafter referred to 
as ADHRS models) 



/ z u C \ 

z' u C y u E x u B 
V x' u B A J 



( z e C \ 

z' e C y e E x e B 
V x' e B A J 



/ z d C C \ 

= z' d y d E x d B , X E 

V x' d B A ) 

(2.2) 

where the x, y, z's are Clebsch factors arising from the VEV's of the adjoint 
fields. Thus realistic fermion masses and mixings can be obtained. 

The radiative corrections to the soft SUSY-breaking operators above Mq 
are now more complicated. From the interactions (2.1) the following soft super- 
symmetry breaking operators are generated: 

XUA)m 2 ki (A)X e3 (A) 4ft, (2.3) 

where 4>i,<fij are scalar components of the superfields, and Xij(A) are adjoint 
dependent couplings, X(A) = X'j^- ■ ■ ■ jf-- After the adjoints take their VEV's, 
the ml e (A) become the usual soft scalar masses. If we ignore the wavefunction 
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renormalization of the adjoint fields (which is valid in the one-loop approxima- 
tion), this is the same as if we had replaced the adjoints by their VEV's all the 
way up to the ultraheavy scale where the ultraheavy fields are integrated out, 
and treated these nonrenormalizable operators as the usual Yukawa interactions 
and scalar mass operators. This is a convenient way of thinking and we will use 
it in the rest of the paper. 

Above the GUT scale, in addition to the Yukawa interactions which give 
the fermion masses 

Q\uU c Hu, QX D D C H D , E C X E LH D , (2.4) 

the operators (2.1) also lead to 

QX qq QH Ua , E c X eu U c H U3) NX nd D c H U3 , 

QX qe LH Ds , U c X ud D c H Da , NX n eLHu, (2.5) 

where H Us , H Da are the triplet partners of the two Higgs doublets H v and H D . 
Each Yukawa matrix has different Clebsch factors associated with its elements, 
so they can not be diagonalized in the same basis. The scalar mass matrices 
receive radiative corrections from Yukawa interactions of both (2.4) and (2.5), 
which, in the one-loop approximation, take the form 

Am 2 Q oc X v X\j + X D X ] D + 2X qq X\ q + X qe X\ e , 
Ami oc 2X\jXu + X\ u X eu + 2X ud Xl d , 
Am D oc 2X d Xd + X nd X nd 

+ 2ALA 

udi 

Am 2 L oc ^e^e + 3AjjA^ + Xl e X ne , 

Am! oc 2X E X% + 3A eu A^ n . (2.6) 

In the minimal 5*0(10) model, scalar mass renormalizations above Mq arise from 
a single matrix X\j. It is therefore possible to choose a "U-basis" in which the 
scalings are purely diagonal. This is clearly not possible in the general models. 
All scalar mass matrices and Yukawa matrices are in general diagonalized in 
different bases. Therefore, flavor mixing matrices should appear in all gaugino 
vertices, including in the up-quark sector (where they are trivial in the minimal 
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models studied in [1-4]). The up-type quark-squark-gaugino flavor mixing is a 
novel feature of the general models. Its consequences will be discussed in Sec. 
IV. Also, the flavor mixing matrices are no longer simply the KM matrix. They 
are model dependent and are different for different types of quarks and charged 
leptons, and are fully described in the next section. 
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III. Flavor Mixing Matrices in General Superymmetric Standard Mod- 
els. 



In this section we set our notation for the gaugino flavor mixing matrices 
in the supersymmetric theory below Mq, taken to have minimal field content. 
We also give general expectations for these matrices in a wide variety of unified 
theories. 

The most general scalar masses are 6x6 matrices for squarks and charged 
sleptons and 3x3 matrix for sneutrinos, 



m 



m 2 UL (Cu + x uH cot P)v v ' 



1 V(Ct + A^cot/3K m 2 



2 / m 2 DL (Cd + A D/ utan/?)v D \ 

m ° ~ Ud + A^tan/^K )' 

2 / m| L (C E + \ E vt&n[3)v D \ 

mE ~ Ud + A^tan/^K m| fl )' 

™l = KJ , (3.1) 



where mfj L , ™?d l -> m u R i m D fl > m E L i m E R are 3x3 soft SUSY-breaking mass ma- 
trices for the left-handed and right-handed squarks and sleptons, and CuiCdi Ce 
are the trilinear soft SUSY-breaking terms. To calculate flavor violating pro- 
cesses, such as fi — > e7, one can diagonalize the mass matrix m 2 E by the 6x6 
unitary rotation matrix Ve and m 2 , by the 3x3 unitary rotation V v , 



m 



E 



V E m 2 E V E , ml = V u m 2 X, (3.2) 



where m 2 E , m 2 are diagonal. The amplitude for fi — > &-/ is given by the diagrams 
in Fig. 1, summing up all the internal scalar mass eigenstates. 

If the entries in the scalar mass matrices are arbitrary, they generally give 
unacceptably large rates for flavor violating processes. From the experimen- 
tal limits one expects that the first two generation scalar masses should be 
approximately degenerate and the chirality-changing mass matrices £ A should 
be approximately proportional to the corresponding Yukawa coupling matrices 
Xa- In this paper we treat the chirality-conserving mass matrices and chirality- 
changing mass matrices separately, i.e., the mass eigenstates are assumed to be 
purely left-handed or right-handed, and the chirality-changing mass terms are 
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treated as a perturbation. This may not be a good approximation for the third 
generation where the Yukawa couplings are large, the correct treatment will be 
used in the numerical studies of Sec. VI. The superpotential contains 

W D Q T \uU c Hu + Q T \ D D C H D + E cT X E LH D , (3.3) 

where Xu, Xd,Xe are the Yukawa coupling matrices which are diagonalized by 
the left and right rotations, 



Xu 


= Vu l XuV Ur , 


x D 


= VS L A D vJ a 


Xe 


= Ve r XeV El . 



The soft SUSY-breaking interactions contain 

QtmjQ + U c] mlU c + D^m 2 D D c + Vm 2 L L + £ ct m|*£ c 

+ Q T CuU c H v + Q T C D D C H D + E cT C E LH D . (3.5) 

Because the trilinear terms should be approximately proportional to the Yukawa 
couplings, we write 

C = Co + AC = AX + AC- (3.6) 
The soft-breaking mass matrices are diagonalized by: 

m£ =U Q m 2 Q *Ul m 2 u = U u m 2 u U{ J , m 2 D = U D m 2 D U}» 

mi =U L m 2 L Ul m|* = U E m 2 *U E , (3.7) 

ACu = V^ACXI, ACd = V D * L AC D V2 R , ACe = V^ACX*- (3-8) 

In the mass eigenstate basis the rotation matrices V, U appear in the gaugino 
couplings, 
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A, = V2g' \ - \eiWl L e L N n {H n 2 + cot w tf n ~ ) + r L W ER e R N n H nS 

7T = 1 

+ ^ cot 9 w V L V L N n H n ~ A 

+ u L Wl L u L N n { l -H n]3 + 1 cot W #™ 3 ) + J L iV n (^ n ~ - ~ cot W #^ 3 ) 

" |^L^^ n # n ~ + -d c L Wl R d R N n H n ~ + h.c. 



3 ^ «AR 11 nB 1 3 

2 

< 1 " 

c=l 

r t JT,/',, _L ,"7 , T/T/t X/'^t 



+ dxW^fiiCxc^cs) + u L W^d L (xlK*~) + h.c] 

+ V2g 3 [u L Wlj L u L g + d L W ] DL d L g + u c L W ] UR u R g + ~cF l WdJr9 + h.c], (3.9) 

whereQ the neutralino and chargino mass eigenstates are related to the gauge 
eigenstates by e.g. B = E n =i H n S N n, w 3 = Eti^^- * + = ELi-^cJSXc 
and 

w** = ulv EL , W Er = U e V Er , W Ul = tf Q V UL , W Dl = U^ q V Dl) 
w Ur = ulv UR , w Dr = U ] d V Dr . 

There are also non-diagonal chirality-changing mass terms 

-C n m d = e T R W* ER (A E + pLtanP)X E Wl L e L v D + e T R V T E AC, E U L e h v D 
+ d[W^ L (A D + /j tan (3) XdW^JrVd + 2£U%AC D U D d&) D 
+ u T L W^ L (A v + n cot P)\uWl R u R Vu + uIU^ACuUuUrVu 
+ h.c. (3.10) 

The lepton flavor violating couplings are summarized in Fig. 2. 

In the rest of this section we discuss the flavor mixing matrices in the min- 
imal supersymmetric standard model, minimal and general SU(5) and SO (10) 
models, with moderate or large tan/3. The results are summarized in Table 1. 

2 Neutrino masses are not discussed here and we choose the neutrino to be in the sneutrino 
mass eigenstate basis. 
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For the minimal supersymmetric standard model, the radiative corrections 
to the soft masses only come from the Yukawa interactions of the MSSM: 

Anig oc XuMi + k^d^di 

Amy oc 2Aj 7 A ; y, 

Am 2 D oc 2X^X0, 

Am 2 L oc X E X E , 

Am| oc2X E X E . (3.11) 

We have assumed a boundary condition on the scalar mass matrices m 2 A oc / 
at Mpl, and 1 represents the possibility that the proportionality constants 
are not universal. For moderate tan/3, X t ^> Xb so that the radiative corrections 
are dominated by \ t . Thus one can neglect the Ad contribution and the only 
nontrivial mixing is Wp L . For large tan/5, X t and are comparable, so will 
lie between XjjXjj and A^A^. Therefore both Wu L and Wo L are non-trivial, 
with KM matrix and one can ignore them. 

For the minimal SU(5) model, there are only two Yukawa matrices, Xu = 
A10, Ad = Xe = A5, and 

Anig oc 3XuX\j + 2/iXDX^, 

Amy oc 3X\jXu + 2kX j ' d Xd, 

Am 2 D (x4X^ D X D , 

Am 2 L (xAX^Xd, 

Am| oc 3XuXlj + 2kX d X ] d . (3.12) 

For moderate tan/3, A t 3> Xb, we have non-trivial mixings for Wu L and W Er , 
as found in @- For large tan/3, A^ can not be ignored, giving non-trivial 
mixings for Wu L and Wjj r - 

For the minimal SO (10) model considered in 0, Q> 

ArriQ oc 5XuXjj + 5kAdA d , 
Amfj oc 5XjjXu + 5kX^ d Xd, 
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Am 2 D oc 5X v Xu + 5kX d X d , 
Am 2 L oc 5Xjj\u + 5k\^ d \d, 

Am| oc 5XuX\j + SkAoA^. (3.13) 

We have non-trivial mixings Wd l , Wd r , We l , and We r for moderate tan/3 and 
non-trivial mixings for all Ws for large tan f3. 

For the general 577(5) or 5*0(10) models, defined in the last section, we get 
non-trivial mixings for all mixing matrices in general. However, in 577(5) models 
with moderate tan/?, the splittings among w? D and m 2 L are too small (because 
they are generated by the small X^(A)) to give significant flavor changing effects. 

One might expect that the mixing in the Wu's are smaller than those in 
the Wd's because of the larger hierarchy in A{/ compared with Xr>. However, 
a given W is the product of a W (which diagonalizes the scalar mass matrix) 
and a V (which diagonalizes the Yukawa matrix). Even if the mixings in VJ/'s 
are smaller than those in Vd's because of the larger hierarchies in Xu, we do not 
have a general argument for the size of mixings in U matrices. This is because 
U diagonalizes (appropriate combinations of) known Yukawa matrices and un- 
known Yukawa matrices appearing above the GUT scale, (2.5). The mixings 
in W and V can add up or cancel each other. Our only general expectation 
is that these new Yukawa matrices have similar hierarchical patterns as Xu or 
Xd- Without a specific model, one can at most say that all non-trivial W's 
are expected to be comparable to Vkm\ the argument that the mixings in Wu's 
should be smaller than is Wp's is not valid. 

In the minimal models at moderate tan f3, the leading contributions to flavor 
changing processes, such as \i — > e^y, involve diagrams with a virtual scalar of the 
third generation. Although such contributions are highly suppressed by mixing 
angles, they dominate because they have large violations of super-GIMjll]]: the 



top Yukawa coupling makes m~ very different from m~, m~. At large tan j3, the 
strange/muon Yukawa couplings get enhanced, so the splitting between m~ and 
m~ increases, leading to potentially competitive contributions to flavor changing 
processes which do not involve the third generation. The importance of these 
new diagrams can be estimated by comparing the contributions to Am^ (in a 
basis where gaugino vertices are diagonal) when the super-GIM cancellation is 
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between scalars of the first two generations (2-1) and third generations (3-1): 



We can see that for large tan/3 (or any tan/5 with small X s coming from the 
mixing of Higgs at Mq i.e., X s (Mg) = ^^^(Mq)), this could be comparable to 
the flavor violating effects from the large splitting of the third generation scalar 
masses. However, for the /i — > in SO (10) models, it does not contribute 
to diagrams which are proportional to m T , (because it does not involve the 
third generation scalars), the dominant contributions are still those diagrams 
considered in [fJJ. For flavor changing processes which do not need chirality 
flipping, such as K — K mixing, and all flavor changing processes in SU(5) 
models, this non- degeneracy between the first two generations is important. 
The above discussion is summarized in Table 1. 




2 



for A 2 = A c , 
for A 2 = A s . 



(3.14) 
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Table 1 







SU(5) 


50(10) 




MSSM 


Minimal 


general 


minimal 


general 




V 




V 




V 


bm\ 


• 


• 


o 


• 


o 


w Ul 


• 


• 


V 


• 


V 


w Dl 






V 


V 


V 


W Ur 




• 


V 


• 


V 


w Dr 






V* 


V 


V 


W El 






V* 


V 


V 








V 


V 





Table 1: Summary table for the flavor mixing matrices: 

8m\ : important effects due to some third generation scalars not degenerate 
with those of first two generations. 

bm\ : non-negligible effects due to non- degeneracy of the scalars; of the 
first two generations. 

Wi : fermion i and scalar i are rotated differently to get to mass basis. 

\J : present for any value of tan/3. 

• : present only for large tan/5. 

o : present for large tan (3, but model dependent for moderate tan (3. 
-: not present. 

* : although present, its effect for moderate tan/5 on flavor violation is 
small due to the small non- degeneracy among different generation scalars. 
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IV. Phenomenology from up-type mixing 



As discussed in the previous section, unlike the minimal models with mod- 
erate tan/3 studied in UJ, ^ [| f| in generic GUT's (for any tan/?) and even for 
minimal GUT's (at large tan/3), we expect mixing matrices in the up sector. 
Having motivated an origin for non-trivial up mixing matrices Wu L(R) 7^ 1, we 
consider some effects they produce. In the following we simply assume some 
W Ul{r) at the weak scale and consider their phenomenological consequences. 
(See however Section V and the appendix for a discussion of the scaling of mix- 
ing matrices from GUT to weak scales.) In particular we discuss D — D mixing, 
corrections to up-type quark masses, contributions to the neutron electric dipole 
moment (e.d.m.) and the possibility of different dominant proton decay modes 
than those expected from minimal models. 
IVa. D — D mixing: 

To get an idea for the contribution of up-type mixing matrices to D — D 
mixing, we follow |12, |H| and employ the mass insertion approximation. The 



bounds obtained from D — D mixing on the 6x6 up-squark mass matrix m 



u 



m 



m 



Ull 
2 

Url 



rn 



m 



u LR 
u RR 



(in the basis where gluino and Yukawa couplings are diagonal) 



are summarized in [13]. For average up-squark mass of m = 1 TeV, they are 



< .04, 



77T 



77T 



(4.1) 



m 



< .06. 



(4.2) 



Consider first (4.1). In the last section we estimated that the contribution to 
m\ 2 from the slight non- degeneracy between the first two generation scalars is 
generically at most comparabale to that from the non-degeneracy between the 
first two and third generation scalars. Thus, for our calculation, we only consider 
the contribution from the splitting between first two and third generation scalars. 
Then, for A = L,R 



m AA\2 



77T 



Wi 



I M3 W U A 32 



2 2 



< 



VA13 W U A32 



(4.3) 



We see that for Ws of the same size as the corresponding KM matrix elements, 
the LHS of (4.1) is of order 4 x 10~ 4 , and the bound is easily satisfied. Turning 
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to (4.2), note that if £ v = AXy, m 



UlR12 



0. However, we expect £ 



v 



A\jj + A£u, with induced in running from Mpl to Mq having primarily a 
third generation component in the gauge eigenstate basis. If all relevant mixing 
matrix elements are of order the KM matrix elements, we expect 



'LRV2 



o 



Amt 



V td V ts 



Again, we see that the bound (4.2) is generically easily satisfied, 



and thus we do not in general expect significant contributions to D — D mixing. 

IVb. Weak-scale corrections to up-type quark masses: 

It is well known that there are important weak-scale radiative corrections 
to the down quark mass matrix proportional to tan/5 |7], || [14], [15|, 16 . In 
general unified models, with non-zero Wjj, there are also importrant weak scale 
corrections to the up quark mass matrix. 

From the diagram in Fig. 3, we have a contribution to up- type masses 
proportional to m t . We find, again assuming degeneracy between the scalars of 
the first two generations, 



Am! 




h{x UL , 3^jj) -|- h( y x UL , S'u^j)] 



Wu L3i W v [h(x tL ,x tR ) - h(x tL ,x v 



(4.4) 



where 



X j 



m i 1 1 \ 
Ml' h ( x > v > 

9 



x-y 



x log x y log y 



1 — x 



i-y 



(4.5) 



L(R)31 



is 



The largest fractional change in the mass occurs for the up quark. If Wi 
comparable to the corresponding KM matrix element, the contribution to 
is not significant. However, if each of the Wu L , R)31 are a factor 3 larger than 
the corresponding KM elements we can get sizable contributions. In Fig. 4, we 

plot ^ - 2* 



in. H^r — — L space, where we have assumed m 



m u R = m u\ m t L 



ni 



tR 



and we have put \W\ 



u L31 \ = \W Ur31 \ = 1/30, (A + ficotj3)/mj= 3. 
Any deviations from these values can simply be multiplied in Am u /m u . In some 
regions of the parameter space it is possible to get the entire up quark mass as 
a radiative effect. 
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IVc. Neutron e.d.m.: 

If we attach a photon in all possible ways to the diagram giving the contri- 
bution to w-quark mass, we get a contribution to the w-quark e.d.m., which is 
proportional to m t for any value of tan f3. Evaluating the diagram, we find 



d u = elFlsin, 



(4.6) 



where 



F=-(^-) m t A + fiCOtP W u Wu W? T 

a 



where 



and 



G2{%t L , Xt R ) — G2(xt L , x UR ) — G2\Xt R i %u L ) + G(x Ul ,x Ur ) 



n i \ 9{x) -g(y) i x 

x-y 



2x- 1 



(4.7) 



[x 2 - 1 - 2xlogx] (4.8) 



Im 



mtW UL31 W* UL33 W Umi W* Umz 



\rntW UL31 W{j m3 W Umi W^\sin 



K (4.9) 



In general we expect a large non-zero sin <p u . If the combination of W s appearing 
in the above is comparable to the combination giving a down quark e.d.m., 
the u-quark contribution will dominate over the <i-quark contribution to the 
neutron e.d.m. considered in by a factor 4m ^ , (the factor 4 comes from 
the quark model result d n = 4/3<id — l/3<i u ). Hence, the neutron e.d.m. may be 
competitive with \x — ► ej and d e as the most promising flavor changing signal 
for supersymmetric unification. 
IVd. Proton decay: 

Finally we turn briefly to the relevance of up-type mixing matrices for 
proton decay; in particular to the important question of the charge of the lepton 
in the final state. We know that upon integrating out the superheavy Higgs 
triplets we can generate the baryon number violating operators ^^(QQ)(QL) 
and j^(EU)(DU) in the superpotential. These operators must subsequently 
be dressed at the weak scale in order to obtain four-fermion operators leading 
to proton decay. The dressing may be done with neutralinos, charginos or 
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gluinos where possible. Since the dressed operator grows with gauge couplings 
and vanishes for vanishing neutralino/chargino/gluino mass, one might naively 
expect gluino dressing to be most important. However, if the up-type mixing 
matrices are trivial, gluino dressed operators can only lead to proton decay with 
a neutrino in the final state. To see this, we examine each operator separately: 
(eu a ) (dbU c )e abc (where a, b, c are color indices) must involve it's from two different 
generations because of the e abc . One of them has to be a u, so the other is a c or a 
t. If there is no up mixing, the up flavor does not change in the dressing process, 
so the final state would have to contain a c or a t. Since m t , m c > m p , this can 
not happen. Next, consider (QQ)(QL) = u a L d b L {u c L eL — d c L VL)e a bc- By exactly the 
same argument as the above, the u a L d b L u c L eL£ a bc operator can not contribute to 
proton decay. Thus, we see that in the absence of mixing in the up sector, gluino 
dressing can only give neutrinos in the final state. However, the above arguments 
break down if up-mixing matrices are non-trivial, since gluino dressed diagrams 
give a significant contribution to the branching ratio for charged lepton modes 
in proton decay. A detailed study of flavor mixing in the up sector (T7J concludes 
that, whether the wino or gluino dressings are dominant, the muon final state 
in proton decay is of greatly enhanced importance. Without the mixings, one 
expects T -£^^k+&) ~ 10~ 3 . The up mixing in general models increases this by 
0(100) making the mode p —>■ K°fi + a favorable one for discovery of proton 
decay. 
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Section V. Large tan/5: Analytic Treatment 

The large tan /3 scenario is interesting for a number of reasons. For moderate 
tan/?, the only way to understand m t 3> mb,m T is to have At ^> A&, A T at the 
weak scale. This gives us little hope of attributing a common origin for third 
generation Yukawa couplings at a higher scale. However, for large tan/? ~ 
the weak scale A t , A&, A T are comparable and the above hope is restored. 



0('- 



(In fact it is realized in SO (10) models like the ADHRS example outlined in 
section VII). For us, this is sufficient motivation to study the large tan/3 case in 
more detail. Also, this case was not studied in || . We shall see that unexpected 
new features arise in the large tan/? limit. 

The largest contribution to the \x — > amplitude comes from the diagram 
with L — R scalar mass insertion (Fig. 5). In the L — R insertion approximation, 
the amplitude for Hl(r) decay is 
a 

Fur) = ^^m T W Erim W Enrr , Wp W* F (A E + utan/3) 

47T COS 2 9 W ^ L ( R h2 ^«(i)3i E HR) 3 3 £ fl(i)33 V r y ' 

G 2 (m£ ,m£ ) - G 2 (m£ , m~ ) - G 2 (m~ , m£ ) + G 2 (m~ , m£ ) 



where 



G 2 (m\, ml) 



G 2 (m\) - G 2 (mp 

2 2 

m\ — m 2 



2 

m 



G 2 (m 2 ) = £(# g + cot 6 w H n ~ s )g 2 — . (5.1) 

n=l \ n/ 

Note, however, that for large tan /3 the L — R insertion approximation may be 
a bad one, since the chirality changing mass for the third generation becomes 
comparable to the chirality conserving masses. A correct treatment will be used 
for the numerical analysis in the next section. We still expect, however, that the 
amplitude to be proportional to We 32 We 31 because of the unitarity of the mixing 
matrices: the sum of contributions from the first two generations is proportional 
to WiiW*j + W 2i W 2j = - W3iW^ for i ^ j, and the contribution from the third 
generation is itself proportional to W 3i W^. 

Two simplifications in the dependence of the \x — > rate on parameter 
space occur for large tan/3. First, since the dominant diagram involves the 
L — R insertion (A + \i tan /3)m t , and since tan/3 is large, the amplitude does 
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not depend on the weak scale parameter A. Second, in the large tan/3 limit, the 
chargino mass matrix is 

X ~\V2M W cos(3 -fi J \ -fi /' ^ ' ' 

and the parameters M 2 ,/a have a direct interpretation as the chargino masses. 
(Note that this assures us that /xtan/3 will likely always be much bigger than 
A] for a tan/3 of 50, the LEP lower bound on chargino mass of 45 GeV tells us 
that /itan/3 > 2 TeV, so for A to be comparable to //tan/3 we must have A > 2 
TeV.) 

In considering \x — > &-) for large tan/3, two factors come immediately to 
mind which tend to (perhaps dangerously) enhance the rate over the case with 
moderate tan f3. 

(i) As we have already mentioned, the dominant contribution to /t — > e-f 
grows with tan/3; the diagram in Fig. 5 is proportional to tan/?, a factor of 900 
in the rate for tan (3 = 60 compared to tan (3 = 2. 

(ii) For large tan/3, A r can be 0(1) and we can not neglect its contribu- 
tion to the running of the slepton mass matrix from Mq to Ms (soft SUSY 
breaking scale). This scaling generally splits the third generation slepton mass 
even further from the first two generations, meaning a less effective super-GIM 
mechanism and a larger amplitude for \i —>■ e^y. 

While both of the above effects certainly exist, there are also two sources 
of suppression of the amplitude for large tan /3, which can together largely com- 
pensate for the above factors: 

(i)' Large tan (3 allows A t to be smaller than for moderate tan/3. There are 
two reasons for this. First, large tan/3 allows vjj to be larger and so At can be 
smaller to reproduce the top mass. Secondly, b — r unification \TB\ is achieved 



with a smaller A T in the large tan/3 regime |7|, ||. Since X t is smaller, a smaller 
non- degeneracy between the third and first two generations is induced in running 
from Mpl to Ma, suppressing the amplitude compared to the moderate tan/3 
case. 

(ii)' In comparing large and moderate tan/3, we must know how the mixing 
matrices Wl,r 3i (appearing at the vertices of the diagrams responsible for /a — > 
&y) compare in these two cases. In the moderate tan/? minimal models discussed 
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in 0, W^R 3i were equal to the corresponding KM matrix elements Vkmzi at 
Mq, and this equality was approximately maintained in running form M G to 
Ms- As discussed in the previous sections, for more general models one expects 
that the Wl(r) m at M G are equal to Vkm3i at M G up to some combination 
of Clebsches. One might then expect (as in the minimal models) that this 
relationship continues to approximately hold at lower scales. In fact for large 
tan (3 this expectation is false. We find that often, the Wi(R)3i decrease from M G 
to Ms, overcompensating for the increased non-degeneracy between the third 
and first two generation slepton masses induced by large A T (point (ii) above). 

In the following, we examine the scaling of these mixing matrices in detail. 
Consider first the lepton sector. The renormalization group equation (RGE) for 
X E (in the following * = j§2£) is 

- dXE - X E [3X E X E + Tr(3X^ D X D + X E X E ) - 3g\ - ~ g 2 } (5.3) 



dt & ^ s--u - * w 5 

giving 

d 18 
- T X E X E = 6X E X E + 2X E X E Tr(3X D .X E . + X E X E ) - (6g 2 2 + - gf)X E X E (5.4) 
dt 5 

-^-X E X E = QX E X E + 2X E X E Tr(3X ] D X D + X E X E ) - (6g 2 2 + ^ gj)X E X E . (5.5) 
dt 5 

These in turn imply that the basis in which A^A^ is diagonal, and the (in general 

different) basis where A^A^ is diagonal, do not change with scale. Consider now 

the evolution of the left handed slepton mass matrix m|. The RGE for m| is 

-^-m 2 L = (m|+2mi )X E X E +2X E m E X E +X E X E m 2 L +2C E C E + Gaugino terms. 
ill 

(5.6) 

In the basis where A^A^ is diagonal, keeping only the A r contribution, the 3i 
entry (i ^ 3) becomes: 

-^™Li = A>Li+2(Ck E ) 3 *- (5.7) 

In this basis, we have m 2 L = W L rri 2 L WL- (Here and in the remainder of this 
section, we abbreviate W El{r) — > W L ( R )). Assuming degeneracy between scalars 
of the first two generations, mf, 3 . = WL3iWl 33 (m 2 L — m 2 h ) = WLSiWl^Am 2 ^. 
Then (5.7) becomes 

-j t {W LM Wl^m 2 L ) = \ 2 T {W L , t Wl,Am 2 L ) + 2(Cks)«- (5-8) 
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For now, we ignore the (CeCe)^ term in (5.8), yielding the solution: 

(W L3i Wl 33 Am 2 L )(M s ) = e- u (W L3l Wl 33 Am 2 )(M G ), (5.9) 



where 



Thus, 



I 



W)- (5-10) 



WlkWUMs) = ^ T ^^j Wl 3i W L33 (M G ). (5.11) 
Similarly, we find 

W mt WUMs) = e ~ 2Ir ^^j W ^ W ks(M G ). (5.12) 

Note that, generically the quantities ^ (Jl) |^°| are smaller than one, since the 
third generation mass gets split even further from the first two generations in 
running from M G to M s . Thus, we find that the W L ( R ) 3i get smaller in magnitude 
as we scale from M G to M s , in contrast with the KM matrix elements V KM 3i, 
which scale as 

Vkm^Ms) = e {h+Ib) V KM ^M G ). (5.13) 

Suppose that at M G the Wl{r) are related to Vkm though some combination of 
Clebsches determined by the physics above the GUT scale. 

Wi {R)33 W L{R)3i (M G ) = z iL{R) V KM 3i(M G ). (5.14) 

This relationship is not maintained at lower scales; instead we have: 

Wl 33 W L3l {M s ) = ^^^^z tL V KM3l (M s ), (5.15) 

WLW^(MS) = ^^e-^^z lR V KM3l (M s ). (5.16) 
l\m R {M S ) 

The dominant contribution to the — > e 7 rate is proportional to 

\Wl 33 W L32 W R33 W mi {M s )\ 2 + \Wl 33 W L3l W R33 W R32 (M s )\ 2 , giving 



Br(ji — > e 7 ) 



Am 2 L (M G ) Am|(M G ) 



2 



e 



x [\z 2l z 1r \ + \z 1l z 2r \ ) 



Am 2 L (M s ) Am R (M s ) 
x Br(n -> e 7 , W L(R)3i Wl {R)33 (M s ) -> IW(M 5 )) 
= e(l ViJ 2 + |^ 2 J 2 ) x 5r(^ -> e 7 , W[ [R)33 W L{R)33 (M s ) -> IW(M S )). 

(5.17) 
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This e represents a possibly significant suppression of the rate for large tan j3. 

At this point, the reader may object: it is true that the Wl(K)3i decrease 
from Mg to Ms, but as already mentioned, the non-degeneracy between the 
third and first two generating is increasing. Which effect wins? We argue that 
in general there is a net suppression. This is easiest to see if in computing the 
/i — > cy amplitude, we use the mass insertion approximation rather than mixing 
matrices at the the vertices (Fig. 6). Although this may be a poor approxi- 
mation, it serves to illustrate our point. (Of course no such approximation is 
made in our numerical work.) From the diagram it is clear that the amplitude 
is proportional to m 2 L32 m 2 R31 (M s ). From (5.7), we see that the rate scales as 

K 32 ™Li) 2 {M S ) = e-^+ 4 '^) (ml 2 ml si ) (M G ), (5.18) 

a net suppression. In the mass insertion approximation, then, the terms ^m^Ms) 
in (5.17) serve to exactly compensate for the increased non-degeneracy between 
m~ L and m~ ; what remains is still a suppression. This, together with (i)' 
invalidates the naive expectation that the theory is ruled out in most regions of 
parameter space due to the enhancing factors (i) and (ii), (although there are 
still stringent constraints on the parameter space). 

The above analysis suggests that individual lepton number conservation is 
an infrared fixed point of the MSSM (whereas individual quark number conser- 
vation is an ultraviolet fixed point). A more complete analysis of scaling for the 
lepton sector and a discussion of scaling in the quark sector is presented in the 
appendix. 
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Section VI. Large tan/3: Numerical Results 

The amplitude for \i — > depends on the 6x6 slepton mass matrix M 2 . 
In the basis where nr^, vo? E are diagonal, we have 

„ fml +D L k \ 

M = L . 6.1 

where in the large tan (3 limit, Di = —(T^ — Qi sin 2 9w)M^ is the D-term contri- 
bution, and kij = fim T tan pWL 3i W R3 j. The amplitude from Fig. 1 for fj, L decay 
is 



where 



G 2 (W hH ), (6.3) 



G 2 {M 2 ) LL G 2 (M 2 ) LR 
G 2 (M 2 ) RL G 2 (M 2 ) RR 
In 0, M 2 was approximately diagonalized by the /xm r tan/3 insertion approx- 
imation, and G2(M 2 ) was calculated using this approximate diagonalization. 
Since here tan/3 is large, we wish to avoid making such an approximation, and 
numerically diagonalize the full 6x6 M 2 . 

Faced with a rather large parameter space, we must decide which param- 
eters to use in our numerical work. We have firstly decided to do our analysis 
only for large tan /3, since the moderate tan /3 scenario has been covered in [0 . 
Secondly, we choose to present our results in a different way than in [[|, where 
the rates for fi — > ej were plotted against a combination of Planck scale and 
weak scale parameters. In our work, we compute fi — > ej entirely in terms 
of weak scale parameters. In particular, we assume that the necessary condi- 
tion for a significant /i — > rate exist at the weak scale, namely non-trivial 
mixing matrix Wl,r u and non-degeneracy between third and first two genera- 
tion slepton masses. In the previous sections, we have shown a possible way in 
which these ingredients may be produced. Our plots for /i — > rates are made 
against low energy parameters, and we separately plot the regions in low energy 
parameter space predicted by our particular scenario for generating \x — > ej. 
This way, our plots are in terms of experimentally accessible quantities and can 
be thought of as constraining the parameter space of the effective 3-2-1 softly 
broken supersymmetric theory resulting from the spontaneous breakdown of a 
GUT. (We use the GUT to relate weak scale gaugino masses.) Our low energy 



23 



plots have no dependence on the physics above the GUT scale, all the model de- 
pendence comes into the predictions for low energy parameters the GUT makes. 
If the predicted region of low energy parameters corresponds to a /i — > rate 
exceeding experimental bounds, the theory is ruled out. 

There is a more practical reason for working directly with low-energy param- 
eters specific to large tan/5: the well known difficulty in achieving electroweak 
symmetry breaking in this regime. Working with high energy parameters, and 
imposing universal scalar masses necessitates a fine-tune to achieve SU (2) x U(l) 
breaking. However, we have nowhere in our analysis made the assumption of 
universal scalar masses, hence the Higgs masses and squark/slepton masses are 
independent in our analysis, and therefore the \i parameter is not tightly con- 
strained by squark/slepton masses. Working with weak scale parameters allows 
us to assume that the desired breaking has occurred without having to know 
the details of the breaking. 

With the aforementioned assumption about the existence of a GUT, and 
assuming degeneracy between the first two generations, the rate for /i — > de- 
pends on the weak scale parameters /i, tan f3, M 2 , m£ , m~ , rr& , m£ , Wi3i, Wr^. 
We know that the amplitude depends on Wl{r)^i simply through the product 
WL 3i Wii 3 j, so for normalization in our plots we put Wl^r^ = Vkmz%- Any devi- 
ation from this can be simply multiplied into the rate. We also fix tan/5 = 60, 
and put m r L(R) = m e L(R) ~ Al(-R)- Next, we use some high energy bias to relate 
m~ L and m~ : we assume that their difference is proportional to M 2 ( as would 
be the case if they started out degenerate and were split only through different 
gauge interactions), so we put m~ L = m~ — rM 2 . In all specific models we have 
looked at, r is small (less than about .2). We find that, as long as r is small, the 
rate has little dependence on its exact value, so we put r = 0, m~ L = m~ = m~. 
We also found that as long as ^ is close to 1, there is little dependence on its 
actual value either, so we put = = A. 

Now, the fi — > e7 rate depends only on /x, M 2 ,rn~ and A, and we have the 
large tan (3 interpretation of fx and M 2 as chargino masses. Fixing m~ = 300 
GeV, we make contour plots of Br(fi — * cy). The rate scales roughly as fn~~ A 
and fi 2 for scalar masses heavy compared with gaugino masses. In Fig. 7, we fix 
\x and plot in M 2 — A space. In Fig. 8, we fix A and plot in // — M 2 space. In 
Fig. 9, we plot the values of A predicted by the GUT against M 2 , for various 
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values of A f (M G ) and A e (Ms) and for two values of b 5 , the gauge beta function 
coefficient above the GUT scale. In Fig. 10, we plot the suppression factor e for 
the same parameter set as in Fig. 9. We see that, over a significant region in 
parameter space, e is small, between 0.2 and 0.01. 

It is clear from Fig. 7 that, with no suppression, a typical value for A of 0.3 
(x300GeV) would give rise to rates above the current bound of Br(fi — > ej) < 
4.9 x 1 1 1 [|l^] - However, from Fig. 10, the suppression from e is seen to be 
typically 20, allowing A's of up to 0.45 (x300GeV). We see that e is crucial in 
giving the GUT more breathing room, as A's of less than 0.45 are more common. 
From Fig. 8 it is also clear that regions of small fi and M 2 (that is, light chargino 
masses) are preferred. Smaller \i is preferred because it decreases the L-R mass 
\im r tan (3, small M 2 is preferred because in the limit that the neutralino mass 
tends to zero, the diagrams Fig. 5 vanish. We also note that smaller /i, M 2 are 
preferred for electroweak symmetry breaking f7|. ||. 

If fi and M 2 are both small, the lightest supersymmetric particle (LSP) 
can be quite light, (but where it has significant higgsino component, it must be 
heavier than 45 GeV in order to be consistent with the precise measurement of 
the Z width), and it annihilates (primarily through its higgsino components) 
through a Z into fermion antifermion pairs much like a heavy neutrino. The 
contribution of the LSP to energy density of the universe Qh 2 then just depends 
on its mass, and the size of its higgsino components, both of which only depend 
on fi and M 2 in the large tan/3 limit. In Fig. 11, we make a plot of Qh 2 in 
/i — M 2 space. We see that it is possible to get fl ~ 0(1) in some region of the 
parameter space. 
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VII. The Example of ADHRS Models 

In this section, we study the ADHRS models || which are known to give 
realistic fermion masses and mixing patterns. These models are specific enough 
for us to do calculations and make some real predictions. Although not neces- 
sarily correct, they are good representatives of general GUT models. We believe 
that by studying them, one can see in detail the general features of generic re- 
alistic GUT models and the differences between them and the minimal SU(5) 
or 5*0(10) models. 

As mentioned in Sec. II, in ADHRS models, the three families of quarks 
and leptons lie in three 16 dimensional representations of SO (10), and the two 
low energy Higgs doublets lie in a single 10 dimensional representation. Only 
the third generation Yukawa couplings come from a renormalizable interaction 



A33I63I63IO. 



(7.1) 



All other small Yukawa couplings come from nonrenormalizable interactions 
after integrating out the heavy fields. These interactions can be written in 
general as 

16^(^)16^10. (7.2) 

The AaS are fields in the adjoint representation of 5*0(10) and their vevs break 
5O(10) down to the Standard Model gauge group. Therefore, these Yukawa 
couplings can take different values for fermions of the same generation with 
different quantum numbers under SU(3) x SU(2) x U(l) and a realistic fermion 
mass pattern and nontrivial KM matrix can be generated. In ADHRS models, 
the minimal number (four) of operators is assumed to generate the up, down- 
type quark and charged lepton Yukawa coupling matrices Xu, \d and A#, and 
they take the form at Mq 

( 
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z e C 
y e E x e B 
x' e B A J 
' (7-3) 



where the x, y, z's are Clebsch factors arising from the VEV's of the adjoint 
Higgs fields A a . This form is known to give the successful relations V u b/V c b = 
^Jm u /m c and V t d/V ts = \Jm d /m s p(J so it is well motivated. Strictly speaking, 
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the interaction (7.2) become the usual Yukawa form only after the adjoints A a 
take their VEV's at the GUT scale. However, as we explained in Sec. II, they 
can be treated as the usual Yukawa interactions up to the ultraheavy scale 
(which we will assume to be Mp^) where the ultraheavy fields are integrated 
out if the wavefunction renormalizations of A a 's are ignored. In the one-loop 
approximation which we use later in calculating radiative corrections from Mp^ 
to Mq, they give the same results, because the wavefunction renormalizations of 
the adjoints A a only contribute at the two-loop order. This makes our analysis 
much easier. Above the GUT scale, in addition to the Yukawa interactions (2.4) 
which give the fermion masses, we have the interactions (2.5) as well. Each 
Yukawa matrix has different Clebsch factors x, y, z associated with its elements. 
All the Yukawa matrices have the ADHRS form 



Ay 



/ ziC \ 

z'jC yiE xiB 
y x'jB A J 



qq, eu, ud, q£, nd, n£. 



(7.4) 



If each entry of the Yukawa matrices is generated dominantly by a single opera- 
tor, like in the ADHRS models, then the phases of the same entries of all Yukawa 
matrices are identical. One can remove all but the A22 phases by rephasing the 
operators. After phase redefinition only E is complex and is responsible for CP 
violation. In order to generate the realistic fermion mass and mixing pattern, 
one expects the following hierarchies, 



V cb ~ e 2 , 



B 
A 
E 
A 
C 

— ~sin# r . ~ e, where e ~ 0.2. 
E 



~ e 



(7.5) 



The hierarchical Yukawa matrices can be diagonalized approximately |j20|| , the 
unitary rotation matrices which diagonalize them at the GUT scale can be ap- 
proximately written as 
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(7.6) 
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The soft SUSY-breaking scalar masses for the three low energy generations and 
trilinear A terms are assumed to be universal^ at Planck scale Mpl as in ||. 
Beneath Mpl, the radiative corrections from the Yukawa couplings destroy the 
universalities and render the mixing matrices non-trivial. In the one-loop ap- 
proximation, the radiative corrections to the soft SUSY-breaking parameters at 
Mq are simply related to the Yukawa coupling matrices and therefore the re- 
lations between general mixing matrix elements and KM matrix elements are 
also simple. This allows us to see the similar hierarchies in the general mixing 
matrices and the KM matrix very clearly. Although the one-loop approxima- 
tion may not be a good approximation for quantities involving third generation 
Yukawa couplings, we will be satisfied with it since it simplifies things a lot and 
the uncertainties in other quantities such as Clebsch factors are probably much 



3 If the nonrenormalizable operators already appear in the superpotential of the underlying 
supergravity theory, the A terms will be different for different dimensional operators, and 
will induce unacceptably large fi — > ej rate because the triscalar interactions and the Yukawa 
interactions can not be diagonalized in the same basis for the first two generations. In theories 
where the nonremormalizable operators come from integrating out heavy fields at Mpl and 
all the relevant interactions have the same A term, the resulting nonrenormalizable operators 
will also have the same A term. 
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bigger than the errors made in the one-loop approximation. The RG equations, 
for m 2 E as an example, from M PL to M G are 



d 9 
—mi 

dt E 



1 



16tt 2 



2(2A £ m|A^ + 2X E m 2 Ho X E + m 2 E X E X E + X E X E m 2 E + 2C E C E ) 



+ 3(2A c „m^Aj u + 2X eu m 2 Hu Xl u + m|A 

euS 



t ) 



gaugino mass contribution 



(7.9) 



In the one-loop approximation, the gaugino mass contributions are diagonal and 
the same for all three generations, so they can be absorbed into the common 
scalar masses and do not affect the diagonalization. The corrections to scalar 
masses at M G have the following leading flavor dependence 



Am| oc 2X E X E + 3X eu X\ u 

( ^c 2 



ZeTeCE* 



\ 



z^CE z' 2 C 2 + y 2 \E\ 2 + x 2 B 2 y e x' e EB + x e BA 

x~?B 2 + A 2 



,(7.10) 



V z e x' e CB y e x' e E*B + x e BA 

where the overline represents the weighted average of the Clebsch factors, z 2 = 
\{2zl + 3z 2 u ) and so on. Because Am 2 E is hierarchical, assuming no big x,y 
Clebsches (ADHRS models have some big z Clebsches), the rotation matrix 
which diagonalizes it can be given approximately as 



Ue(M g ) 



-S El e~^ 



\ S El S E2 e~ 



'£3 



s El 

-S E2 e 



Se3 ^ 

S E 3 

1 / 



(7.11) 



where 



Se 2 — 



x P B 



, s 



z e xLC B 
A 2 



z e y e C\E 


„-i4> _ E 


z' 2 C 2 + y 2 


E 


2 + {x 2 e -x 2 \B 2 ' " \E\ 



Similarly, for other scalar masses the leading flavor dependent corrections at M G 
are 



Am£ oc X E X E + 3A^A^ + X nl X n ^ 
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Am 2 Q oc \uX\j + X D X D + 2\ qq \l q + X qe X\ e , 
Ami oc 2XIXU + X\ u X eu + 2X ud Xl d , 

Aui 2 D oc2X j D X D + Xl d X nd + 2Xl d X ud , (7.12) 

and the rotation matrices which diagonalize them are given by expressions simi- 
lar to (7.11) with Clebsches replaced by the appropriate ones. Then, the mixing 
matrices appearing at the lepton-slepton-gaugino vertices are given by 



w El (m g ) = uIv El 

( 1 _ S EL1 -S L1 e^-^ S Ll (S L2 -S EL2 )e^-S L3 \ 

(S Ll - S ELi )e-^-^ e^-M -(S L2 - S EL y* 

\S ELi (S L2 -S EL2 )e-^ + S L3 (S L2 -S EL2 )e-^ 1 J 

(7.13) 

W El (M g ) = U e V Er 

( e**° _ S Em e^e-S El e^ S El (S E2 - S ER2 )e^ - S Es \ 

~ S El e^ - S Em e^ e* -(S E2 - S^e'^ 

X—SEmiS^—Ss^ + Sfye 1 ^ S E2 — S Er2 1 / 

(7.14) 

where 

C Z (P Q Ze ^ rpi L, 771 X e X e B 2 

Ll = ~E r ' Rl = _ E r ' e = l2/ a — I' 



X ' P B „ x e B , / . ^. x'x e B 2 ^ 

S EL2 = S ER2 = -Sp e = art? U |£| e * - 

- (if J/e ~a; e ,0, 



5 _ feyegjigj = _ Xj3_ - Z e x' e CB 

El ~^C* + f e \E\ 2 + (£* - x 2 e )B* ' E2 = A ' ^ " A* ' 

Ll ^ + ^ + (3-^'^ A ' i3 ~ ^2 > 
and 

x e = — (2x e + 3x eu ), x' e = ^(^e 4" 3x^£ + x n f) etc. 
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Note that 

S L3 ~ Clebsch x S ELi (S L2 - S EL2 ) = Clebsch x (7.15) 

If there is no very big or small Clebsch involved and no accidental cancellation, 
Sl 3 ,Se 3 can be neglected in W's. 
Compared with Vkm, 

V K m(M g ) = vIv Dl 



i(4>d-4>u) 
\S Dli (S Dli - S UL2 )e-»< 



Su Ll - Sd Li z 



e -i(<Pd-<Pu) 

-{S D , -S UL )e 



-S Uli (S DLr 
(S DL2 -S UL2 )e^ 
1 



(7.16) 



where 



s u L1 
S Dli 

SDt.2 



z C 
z d C 

x u B 
A ' 

x d B 
A ' 



VuE - 
VdE- 



Xy,X u B 



A 

x d x' d B 2 



I / 771 x u x u B 

, <Pu = arg I y u E — 

x d x' d B 2 ^ 



b d = arg y d E 



A 



we can see that the IV's and Vkm do have similar hierarchical patterns, but 
have different Clebsch factors associated with their entries. 

When a specific model is given, one can calculate all the Clebsch factors 
and make some definite predictions for that particular model. For example, the 
ADHRS Model 6, which gives results in good agreement with the experimental 
data, has the following four effective fermion mass operators 

33 = 16 3 10 16 3 , 



O 



23 



Ay Ay 

16 2 10 -fl6 3 , 



23 = 16 2 -r-r 



162^- 10 —7—^ I62 or other 5 choices, 
M Ax 



M ) 



M ) 



(7.17) 
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where A X ,A Y ,A B _ L are adjoint's of SO (10) with VEV's in the SU(5) singlet, 
hypercharge, and B — L directions. There are six choices of O22 operators 
which give the same predictions for the fermion masses and mixings, but dif- 
ferent Clebsches for other operators appearing above M G . Fortunately, they 
do not enter the leading terms of the most important mixing matrix elements 
We L32 ,We L31 ,We R32 ,We R31 , which appear in the leading contributions to the 
amplitudes of LFV processes and the electric dipole moment. 

The magnitude of the mixing matrix elements Vkm32, Vkm3i, W Bl32 , We L31 , 
W Br32 , We R31 , and the relevant Clebsch factors are listed in Tables 2 and 3. 
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Table 2: Clebsch factors for Yukawa coupling matrices in ADHRS model 6. 
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Table 3 





ADHRS models 


Model 6 


Relevant process 
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ZdVe{x d -X u ) 2 










z' e yd(x e -x e )(x' e -x' e ) 
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V2z' e z e y d 
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d e 


Wl^3 2 ^ 31 l 2 +I^H32^3ll 2 / Vl4j 


yj (z 2 +z' 2 )z d y e 



Table 3: The relevant Clebsch factors for /x — > ey and c? e in ADHRS 
model 6. 



In ADHRS models tan /3 is large. The fi —>■ ej rate for large tan j3 has been 
calculated in Sec. V and VI for We L32 = We R32 = V ts and We L31 = We R31 = Vtd- 
To obtain the predictions of ADHRS models we only have to multiply the results 
by the suitable Clebsch factors. The relevant Clebsch factors for Model 6 are 
listed in Table 3. For a generic realistic GUT model with small tan/?, for example 
the modified ADHRS models in which the down type Higgs lies predominantly 
in some fields which do not interact with the three low energy generations and 
contain only a small fraction of the doublets in the 10 which interact with the 
low energy generations [21], most of analysis should still hold. In this case the 
leading contributions to fi — > ej are the same ones as in the minimal SO(10) 
model of Ref. Q (Fig. 10 b L)R , c L<R , c' L R of gj). The diagrams c LR ,c' LR involve 
the corrections to the trilinear scalar couplings. 
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In the one-loop approximation the leading corrections to £ E at Mq contain 
pieces proportional to X E , Xe(X e \ e + 3A^A^ + A^A^), (2A E A^ + 3A eu A* JA E 
respectively. The piece proportional to Ag can be absorbed into £ Eq by a redef- 
inition of A E , the other two pieces are proportional to the product of X E and 
the corrections to the scalar masses, 

AC E = AC Efl + AC Ei 

&Ce r = —Am 2 E \ E 

Ve r 

AC El = —\ E Aml (7.18) 

where \i El1 [i Er are proportional constants (\x Er = ^e l = 6m °^ A ° in one-loop 
approximation). The LFV couplings in Fig. 2e, e^U E A^ E U E eLVD, now can be 
written as 

-^—e^U E Am E X E U L e L v D + ——e\U E X E Axv? L U L e L v D 

V-Er He l 

1 _rn . n + 1 



VEr " ~ ~" ~" HE L 



e L R Am E W ER \ E Wl h e L v D + e^W ER \ E Wl L Ara.ie L v D , (7.19) 



where the overline means that the matrix is diagonal. Again, the amplitudes are 
given by the same formulas as in (eqn. 29, 30), except that V 3 %V 3 \(Vi 3 *) 2 has 
to be replaced by W El32 W Er31 W El33 W Er33 , and W Er32 W El31 W Er33 W El33 , and 
|/q by —^-Am E33 and j^-Afn L33 . The results in are only modified by some 
multiplicative factors and therefore represent the central values for the LFV 
processes. 

It was pointed out in 0, |3| that the electric dipole moment of the electron 
(4) constitutes an independent and equally important signature for the SO (10) 
unified theory as \x — > ey does. The diagrams which contribute to the electric 
dipole moment of the electron are the same as the ones which contribute to 
/i — > ej, with \i E (\x c L ) replaced by e E (e c L ). Thus a simple relation between d e and 
the \i — > e7 rate was obtained in the minimal 5*0(10) model [0, 

a 

—it-,, i 
2 M ' 



T(fi^e 1 ) = -ml\F 2 \ 2 , (7.20) 
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td 
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sin0, (7.21] 
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where is an unknown new CP violating phase defined by 

Im{m T (V 3 \nV£) 2 } = |m r (y 3 e i) 2 ra 2 |sin0. 

In a more generic 5*0(10) model, such as the ADHRS model, we still have this 
simple relation but the mixing matrix elements have to be replaced by the W's: 



2r(/^e 7 ) V2\W Emi W, 



E 



^ 1 am l \I\We L32 W Er .^ + \W Er31 W^ 
where 0' is defined by 



R3i I ™„ aJ 



sin0', (7.22) 



Im[m T W EL31 W Emi W Em: W EL J = \m T W Et31 W ERZ1 W Em3 W Em3 \ sin0'. 

In particular, in ADHRS models there is only one CP violating phase, so the 
phase 0' can be related to the phase appeared in the KM matrix of the Standard 
Model. From eqn. (7.13) (7.14) (7.16) we can see that 0' ~ e , e ~ 0^ ~ 0, n = 
(because y u = 0). The rephrase invariant quantity J of the KM matrix is given 
by 

j = imv ud v tb v; d v: b 

~ -S Uli S Dli {S Dl2 - S UL2 ) 2 sm<j) d . (7.23) 

Therefore the CP violating phase appeared in d e related to the CP violation in 
the Standard Model by 

sin0'^ J (7.24) 

\Vtd\\Vub\ 

Finally, as mentioned in the Sec. Ill, we consider the possibility that the 
slight non- degeneracy between the first two generation scalar masses could give 
a significant contribution to the flavor changing processes because of the larger 
mixing matrix elements. We still use ADHRS models as an example to estimate 
this contribution to the LFV process /i — > cy. For an order of magnitude 
estimate, the mass insertion approximation in the super-KM basis employed in 
]1|] will serve as a convenient method. After rotating the Am 2 E in eqn. (7.10) 
to the charged lepton mass eigenstate basis, the contribution from the first two 
generations to Am E2i is 

Am| 2i (2-1) ~ V ER22 V Em Am E22 + V Em2 V Em Am En + V ER22 V* m Am 2 E21 
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z e C z' 2 C 2 + y*\E\* + x 2 B* z e Cz 2 C* i(f> 

"r 1-1. i ^ ZeVt 



E' 



A 2 



E' A 2 



CE, . o 



y2 C \E\+X^ + Z-r e C\E\ 

A 2 



(assume z e = z' e as in ADHRS model). (7.25) 
Compared with the result found in [|l| for minimal SU(5): 
AmljBH) =V*V td Am 2 E33 



ZdC (x d 



2 B 2 



E> 



A 2 



Z d { x d — Xu) 2 ^§\- 



Vd 



Am E 33 



A 2 



(7.26) 



we can see that if the Clebsch factors are 0(1), this contribution is comparable 
to that of the minimal SU(5) model. In order for this contribution to be com- 
petitive with the dominant diagrams (Fig. 10 bL,R,CL,R,c' LR of |2[]) which are 
enhanced by large Clebsch factors are required. While it is possible to have 
large Clebsch factors, we consider them as model dependent, not generic to all 
realistic unified theories. 
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VIII. Conclusions 

In supersymmetric theories, the Yukawa interactions which violate flavor 
symmetries not only generate the quark and lepton mass matrices, but neces- 
sarily also lead to radiative breaking of flavor symmetries in the squark and 
slepton mass matrices, leading to a variety of flavor signals. While such effects 
have been well studied in the MSSM and, more recently, in minimal unified 
models, the purpose of this paper has been to explore these phenomena in a 
wide class of grand unified models which have realistic fermion masses. 

We have argued that, if the hardness scale A H is above M G , the expec- 
tation for all realistic grand unified supersymmetric models is that non-trivial 
flavor mixing matrices should occur at all neutral gaugino vertices. These addi- 
tional, weak scale, flavor violations are expected to have a form similar to the 
Kobayashi-Maskawa matrix. However, the precise values of the matrix elements 
are model dependent and have renormalization group scalings which differ from 
those of the Kobayashi-Maskawa matrix elements. 

It is the non-triviality of the flavor mixing matrices of neutral gaugino 
couplings in the up quark sector which strongly distinguishes between the general 
and minimal unified models, as shown in Table 1. Although the minimal unified 
models provide a simple approximation to flavor physics, they are not realistic, so 
we stress the important new result that flavor mixing in the up sector couplings 
of neutral gauginos is a necessity in unified models, this leads to four important 
phenomenological consequences. While the D° — D° mixing induced by this new 
flavor mixing is generally not close to the present experimental limit, it could 
be much larger than that predicted in the standard model. 

The new mixing in the up-quark sector implies that there may be significant 
radiative contributions to the up quark mass matrix which arise when the super- 
partners are integrated out of the theory. This is illustrated in Figure 4, where 
the new mixing matrix elements have been taken to be a factor of three larger 
than the corresponding Kobayashi-Maskawa matrix elements. In this case the 
entire up quark mass could be generated by such a radiative mechanism: above 
the weak scale the violation of up quark flavor symmetries lies in the squark 
mass matrix. 

The electric dipole moment of the neutron, d n , is a powerful probe of 
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the neutral gaugino flavor mixing induced by unified theories. In the mini- 
mal SO (10) theory, d n arises from the flavor mixing in the down sector, which 
leads to a down quark dipole moment, dj. However, in realistic models the 
flavor mixing in the up quark sector leads to a d u which typically provides the 
dominant contribution to d n . Thus the neutron electric dipole moment is a more 
powerful probe of unified supersymmetric theories than previously realized. 

The presence of flavor mixing in the up sector plays a very important role 
in determining the branching ratio for a proton to decay to K°fi + . In the 
minimal models, without such mixings, this branching ratio is expected to be 
about 1CT 3 : the charged lepton mode will not be seen and experimental efforts 
must concentrate on the mode containing a neutrino, K + u. However, including 
these mixings the charged lepton branching ratio is greatly increased to about 
0.1. While this number is very model dependent, we nevertheless think that this 
effect greatly changes the importance of searching for the charged lepton mode. 

These four phenomenologocal consequences are sufficiently interesting that 
we stress once more that they appear as a necessity in a wide class of unified the- 
ories. The absence of mixing in the up sector is a special feature of the minimal 
models. Since the flavor sectors of the minimal models must be augmented to 
obtain realistic fermion masses, any conclusions based on the absence of flavor 
mixings in the up sector are specious. 

A second topic addressed in this paper is the effect of large tan (3 on the 
lepton process, /x — * which is expected in unified supersymmetric SO (10) 
models. The amplitude for this process has a contribution proportional to tan (5. 
In this paper, we have found that the naive expectation that large tan f3 in su- 
persymmetric SO (10) is excluded by /i — > is incorrect, at least for all values 
of the superpartner masses of interest. Contour plots for the /j, — > branching 
ratio are shown in Figures 7 and 8. It depends sensitively on the parameter 
A, which is the mass splitting between the scalar electron and scalar tau, and 
is plotted in Figure 9. Lower values of the top quark Yukawa coupling, which 
for large tan/3 still give allowed predictions for the b/r mass ratio, give a much 
reduced value for A, thereby reducing the \x — > rate and partially compen- 
sating the tan 2 /? enhancement. A further significant suppression of an order of 
magnitude is induced by the renormalization group scaling of the leptonic flavor 
mixing angles, and is shown in Figure 10. The net effect is that while the case 
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of tan/3 pa m t /m h is not excluded in SO(10), the /x — > e7 rate is still typically 
larger than for moderate tan/3, so that this process provides a more powerful 
probe of the theory as tan (3 increases. 

For large tan (3, \x and M 2 become the physical masses of the two charginos. 
The /x — > e7 contours of Figure 8 show that /i and M 2 should not be too large, 
providing an important limit to the chargino masses in the large tan/3 limit. 
Furthermore, this constrains the LSP mass to be quite small. We find that 
in this region it is still possible for the LSP to account for the observed dark 
matter, and even to critically close the universe, as can be seen from Figure 11. 
However, the requirement that the LSP mass be larger than 45 GeV suggests 
that the two light charginos will not be light enough to be discovered at LEP II. 

As an example of theories with both a realistic flavor sector and large tan (3 
we studied the models introduced by Anderson et al. The flavor sectors of these 
theories are economical: the free parameters can all be fixed from the known 
quark and lepton masses and mixings. Hence the flavor mixing matrices at all 
neutral gaugino vertices can be calculated. These are shown for the lepton sector 
of model 6 in Table 3. The Clebsch factors enhance the /x — > cy amplitude by 
a factor of 2.3, and suppress d e by a factor of 3. Even taking the top quark 
Yukawa coupling to have its lowest value the rate for /x — > in this theory 
is very large. Another interesting feature of these theories is that the flavor 
sectors contain just a single CP violating phase. This means that the phase 
which appears in the result for d n and d e can be computed: since it is closely 
related to the phase of the Kobayashi-Maskawa matrix it is not very small. That 
which appears in d e is given in eqn. (7.24) and is numerically about 0.2. We 
have computed the radiative corrections to m u in the ADHRS models and have 
found that the new mixing matrices in the up sector are not large enough to yield 
sizable contributions: thus the ADHRS analysis of the quark mass matrices is 
not modified. Furthermore, due to a cancellation special to these theories, there 
is no contribution to d n from the up quark at one loop. 
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Note Added: 

While finalizing this work, we received a preprint by Ciafaloni, Romanio 



and Strumia p2| , where the large tan/5 scenario is also considered. However, 
unlike this work, they assume strict universality in soft scalar masses, such that 
imposing electroweak symmetry breaking leads them into a region of parameter 
space with a high mass (1 TeV) for the sleptons. In their discussion of general 
models, they do not include flavor violating RG scaling of scalar masses above 
M G . 
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Appendix 

In this appendix, we first give a more complete treatment of mixing matrix 
scaling in the lepton sector, and then give a treatment for the quark sector. 

Let us return to (5.7) and consider the effect of including the (C E C E )3i term. 
In general the scaling from M PL to M G will generate a CeCe n °t diagonal in 
the same basis as A^A^, so we expect some non-zero (CeCe)^- From the RGE 
for £ E , neglecting gauge couplings, 

-j t C E = C s [5AtAB+rr(3AU^+AtA E )]+A B [4AtC£;+rr(6C J3 Aj 3 +2C £; At)]. 

(Al) 

We have 

-j t (C E C E ) = 5[CUe^e + Xe^eCeCe] + 2Tr(3X^ D X D + X e X e )CUe 

+ 8C e X e XUe + (CLAs + AUs)Tr(6C fl A t D + 2C E X E ). (A2) 

Then, to first order in the off diagonal parts of CeCe an d CeGei anc ^ keeping 
only third generation Yukawa couplings we have 

-j t (CU E h = (CUeU^XI + 6\t + 6^A 6 A T ], (A3) 

where 77 = ^ S3 -- Because of the large numerical coefficient in front of A^, Aj; in 
the above equation, (CLCb)3« is driven to zero more rapidly than Wlm, after 
which it ceases to have any effect on the running of Wlm- More explicitly, from 
(5.7) we have that 

jyL*{t)eti«*W) = -2(CU E h(t)ef>'^'\ (A4) 

Solving (A. 3) for (CeCe)^) an d inserting into (A. 4) we get 



j t (ml 3l (t)ei:> dt ' x ^) = 2(cLC E )3.(M G )e-/o^[^ + 6A ?+ 6,A 6 A t ]( t '). (A5) 



_ d 

Integrating (A. 5), we find 

m 2 L3i (M s )e^ - mi 3t (M G ) = -2 dt& - Jl^x^l+^W) 

Jo 



x (CUeUMg) 

= 5(CUeUM g ). (A6) 
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So, we have 

m 2 L3i (M s ) = e-^[m 2 L3i (M G ) + 5(CUe)(M g )}. (A.7) 

We expect m 2 L3i and (CeCe)^ to be related by some combination of Clebsches 
x at Mq as follows: 

(CUs)* = hxm 2 LU (A8) 

Where A ,ml are the universal A parameter and scalar mass at Mpl, respec- 
tively. Then, we have from (A.7) 

Wt 33 W L3t (M s ) = e-^ ^/fFh l + S^ 2 x]Wt 33 W L3l (M G ). (A.9) 

Am z (Ms) m,Q 

Clearly if 5— %x <C 1, inclusion of the (C E C e)si term in (5.7) do not change any 

A 2 

of our results. If 5^x ~ 1 or > 1, we can still of course use (A.9), but the 

suppression effect may disappear. A simple estimate shows, however, that 5 

i A 2 
itself is already small ~ ^, and so we are only in trouble if -^x is big. To see 

this, replace A T , \b and r\ by some average values A T , \ and fj in the expression 

(A.6) for 5. Then, 



Mq 

S = -2 f™* ° STT * e -t(16\ 2 T +6\ 2 +6\ b \tv) 

Jo 



8\ 2 + 3(\ 2 + fj\ b Xr) 



^log^(16A?+6A2+6^ 6 A r ) 



(A10) 



So, 



151 < 8A?T3(AiTTO (AU) 
For the A's between 0.5 and 1, and fj ~ 1, \5\ ranges from | to ^. 

How can we qualitatively understand the above results for the scaling of 
mixing matrices? The renormalization group equations try to align the soft 
supersymmetry breaking flavor matrices with whatever combination of flavor 
matrices responsible for their renormalization. However, because a given cou- 
pling can only be renormalized by harder couplings, there is a hierarchy in which 
flavor matrices affect the running of others. The Yukawa matrices, being dimen- 
sionless, can only be affected by other Yukawa matrices. In the lepton sector, 
this is the reason that the basis in which e.g. A^A# is diagonal does not change. 
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Next, the soft trilinear terms, having mass dimension one, can only be affected 
by other trilinear terms and Yukawa couplings. Again in the lepton sector this 
means that e.g. CeCe tries to align itself with \%\e- Finally, the scalar mass, 
having dimension two, are affected by everything: m L tries to align with A^A^, 
but suffers interference from CeCe-: unless CeCe * s diagonal in the same basis as 
A^A^. Even if CeCe * s n °t diagonal in the same basis as A^A^, it is trying to 
align itself with A^A S , so m| will still tend to align with A^A^. 

From the above discussion, it is clear that the situation is slightly compli- 
cated in the quark sector. In the lepton sector, there was a fixed direction in 
flavor space given by Xe, with which the soft matrices aligned. In the quark 
sector, we have both A^ and Xp, and A^A^, A^A^, are misaligned (Vkm !)■ 
This complicates the analysis for Wu l ,Wd l so we discuss them last. Let us 
now examine the scaling of W Ur , W Dr . (Throughout the following, we assume 
degeneracy between first two generation scalar masses, we neglect all Yukawa 
coupling matrix eigenvalues except those of the third generation, and we do not 
include the effect of trilinear soft terms in the scaling. The last assumption is 
made for simplicity; we can make similar arguments about the importance of 
these neglected trilinear terms as we did above in the lepton sector.) 

First, we show that the basis in which X\jXu is diagonal remains fixed. The 
RGE for \\j\u is 

-j^Uu = 6(A[ / A^) 2 + 2A[ / A D Al ? A^ + 2(3TrA^A[ / -^3 2 -3^-^ 1 2 )A[ r A C7 . 

(A12) 

Working in a basis where is diagonal, let us see if ^X\jXu has off-diagonal 

components. We have, (recalling that in this basis X^ d Xd = VkmX 2 d Vk M ), 

— ~n(^u^u) a — ^(X-uVkmXj^V^Xu)^ 

at i^j 

= 2X u y K M ll x^ De v^ Ml] X Uj 

= 0fori,jV3 (A 13) 

since we neglect all Yukawa's except the third generation. Similarly, the basis in 
which A^Ao is diagonal does not change. Thus, the discussion for the scaling 
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of W Ur , W Dr is completely analogous to that in the lepton sector, and we find 



W DR3i Wh R jM s ) = e ~ 2Ib ^^W Dm Wh R3 (M G ). 



(A.U) 



(A15) 



We now turn to Wu L , We l - Let Vj} L (t) be the matrix diagonalizing A^A^t): 

\ u \Ut) = V^ L (t)X 2 u (t)V*l(t). (A16) 

In the superfield basis in which XuXjj is diagonal, the squark mass matrix is 
™ 2 Q3i = Vu L m Q V u L - Note as before that m 2 Q * 3i = (Wlm 2 Q *W UL ) 3t = W Ul3 WI lx A 
so we are interested in J^m^. Now, 



dt 



dt 



d 

dt Xu > 



d 



d 



dt 
(A17) 



The second term is the analogue of what we have already seen in the lepton and 
right-handed quark sector; using the RGE for nig we find to leading order 

(<l m Q^) 3 . = ~^ + A ') ffi Q3,- (A18) 
Now, Vu l -^V Ul is obtained from the RGE for A^Ajy. Actually, note that 



< 3>Aj, V Ul 



dt 



» 2 



—> 2 

dt 



+ -y v , (A19) 



so that only \Yu L ^Vu L , ^u] is determined. (This is a reflection of the fact that 
Vu L is not unique: let X(t) be any unitary transformation leaving Mg(t) in- 
variant: m^i) = X\t)rn 2 Q{t)X{t). In our case, X(t) is most generally a U(2) 
matrix in the first two generation subspace. Then, if Vjj l diagonalizes nig , so 
does V Ul X. Under this change, ^Vu L is not invariant, but [Vu L jjVu L , 
is invariant). Further, since we neglect first two generations Yukawa eigen- 
values, [Vi L f t V UL , = for i,j = 1,2, and only [V^ L f t V UL , Aj] 3i(i3) = 
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(t)^'^u l -^V^ is determined, and we can choose all other components of 
V u L ft V u L to vanish. From the RGE for \uX\j, 

-j t {\u\ ] u) = H^uXh) 2 + 2(3TrXu\l - fg 2 3 - 3g 2 2 - ^Xu^u 

+ {X u Xl,X D X t D }, (A20) 



we find 



and thus 
Thus to leading order 



- {^u L {-^uX\j)Vu^ — {X 2 V , V K MX 2 D V^ M } 3i 

= X 2 \lV K M33V KM3i , (-4-21) 

[ V Lj t V ^ = ->^M M3i V KM33 . (A.22) 



[Vi L m 2 Q V UL X L jVu L h = -Am^ A, 2 V^ M3i V KM33 , (A23) 
and finally we have 

-J t (Wu L3 Wl L3 Aml) = (A 2 + Xl)W UL3 Wl L3 A m2 Q + XM M3l V K M3 3 Am 2 Q . 

(A24) 

Similarly we find 

-J t (W DL3i Wl L33 ^m 2 Q ) = (A 2 + \ 2 b )W DL3i Wl L3 Am 2 Q + \M M3l V KM33 Am 2 Q . 

(A25) 

We can formally solve the above equations, e.g. 

-(^+fM Gd ^l V " M3lV T m ) A W 2 (M \ 
W UL M L33 (Ms) = e V ) ^^W Ul3 WUMo), 

(A26) 

and, to a good approximation, given that Wu L3i does not scale very significantly, 
we can replace 

[ MG dt'X 2 b V ^ VK r 3 « h V « M3 * VK r 3 (M G ). (A27) 

Jm s b w UL . M wl L33 w UL3 wl L33 
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So, an approximate solution of the RGE for W Ul , W Dl is 



(A28) 
and similarly 



/ f , + M L+ „ f> - u: "„',V"- (-W ( .-) I I Am 2 

Am^(M s ) 



* ^ w* - VJ " G ' j I Am 2 (M G ) + , N 



(A.29) 

The above results are in agreement with qualitative expectations; the extra 
terms in the exponential of (A.28) and (A.29) are a reflection of the fact that the 
bases in which Xu^u and Xd^d are diagonal change with scale. For moderate 
tan P, however, we expect that the basis in which Af/Ajy is diagonal should not 
change with scale, and in this limit the extra term drops out of (A.28). 
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Figure Captions 



Fig. 1. Feynman diagrams contributing to /i — > cy. 

Fig. 2 Lepton flavor violating couplings in general supersymmetric Standard Mod- 
els. 

Fig. 3 Corrections to the up-type quark mass matrix, proportional to m t . 

Fig. 4 Contours for in ^ — — plane, assuming m~ = m~ = m~, mr = 

m~ tR m~ v W UL31 = W UR31 = i, = 3. 

Fig. 5 The diagram which gives the dominant contribution to ji — > in the large 
tan/5 limit. A photon is understood to be attached to the diagram in all 
possible ways. 

Fig. 6 The dominant diagram (for fi — > e^y) in the mass insertion approximation. 

Fig. 7 Contours for Br(fi — > cy) in M 2 — A plane with m~ = 300GeV, A = 
m e i(fl) " mr iW , W BiWM = 0.04, W El(r)31 = 0.01, for (a) » = 100 GeV, 
(b) \x = 300 GeV. Contours for negative \x are virtually identical. To get 
Br\\i — * ej) prediction from a GUT, multiply by appropriate Clebsch, 
and e factor (Fig. 10). 

Fig. 8 Contours for Br(/i -> e7) in ^ - M 2 plane for (a) A = 0.25, (b) A = 0.5, 
with other parameters same as in Fig. 7. The blacked out regions are ruled 
out by the LEP bound of 45 GeV on chargino masses. The thick dashed 
lines are contours for a 45 GeV LSP mass. 

Fig. 9 Plots of the averaged difference between the third and the first two gen- 
erations charged slepton masses A = Al 1 2 Ah , A-l(r) = m 7 L(R) ( at Ms), 
against M 2 , for \{jr& + m 2 ~ ) = (300 GeV) 2 , A t = X b = A T (at M G ) = 
0.5, 0.8, 1.1, A e (M s ) = 1,0,-1, two values of the gauge beta function 
coefficient 6 5 between M G and M PL , (a) 6 5 = 3 (asymptotically free), (b) 
6 5 = —20. Scalar masses are assumed degenerate at M PL = 2.4 x 10 18 
GeV. Mq is taken to be 2.7 xlO 16 GeV. 

Fig 10 Plots of the suppression factor e against M 2 , with the same parameters as 
in Fig. 9. 
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11 Contours for Vlh 2 in fx — M 2 plane in the large tan (3 limit. Dashed lines are 
LSP mass contours of 30, 45, and 60 GeV. For all regions of ttilsp < 45 
GeV in this plot, the Higgsino components of LSP are too big and therefore 
they are ruled out by the Z width. 
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